Local waiting time fluctuations along a randomly pinned crack front 
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The propagation of an interfacial crack along a heterogeneous weak plane of a transparent Plexiglas 
block is followed using a high resolution fast camera. We show that the fracture front dynamics 
is governed by local and irregular avalanches with very large size and velocity fluctuations. We 
characterize the intermittent dynamics observed, i.e. the local pinnings and depinnings of the crack 
front which trigger a rich burst activity, by measuring the local waiting time fluctuations along 
the crack front during its propagation. The local front line velocity distribution deduced from the 
waiting time analysis exhibits a power law behavior, P(v) oc v~ v with n = 2.55±0.15, for velocities v 
larger than the average front speed (v). The burst size distribution is also a power law, P(S) cx 
with 7 = 1.7±0.1. Above a characteristic length scale of disorder Ld ~ 15pm, the avalanche clusters 
become anisotropic, and the scaling of the anisotropy ratio provides an estimate of the roughness 
exponent of the crack front line, H = 0.66, in close agreement with previous independent estimates. 

PACS numbers: Pacs numbers: 62.20.Mk, 46.30Nz, 61.43.-j, 81.40.Np 



The physics community has recently paid a lot of at- 
tention to the study of damaging processes Q, 0,01 • This 
interest is motivated not only by the practical benefits to 
many engineering domains, but also from a more funda- 
mental point of view, by the diverse challenging questions 
brought forward, in particular in statistical physics 
The role of heterogeneities during crack propagation is 
of central importance since they induce local pinnings 
of the crack front and subsequently trigger a very com- 
plex history of the fracture in the material. One of the 
consequences of this phenomenology is the roughness of 
fracture surfaces left by the crack. Indeed, cracks in het- 
erogeneous media exhibit a self-affine morphology, with 
long range correlations. The associated roughness expo- 
nent was found to be very robust for different materials, 
over a broad range of length scales [E IS III li|Jn| ITTj , 
and was further conjectured to be universal 0, la- A 
recent work 0, suggests that the origin of these self- 
affine long range correlations comes from the elastic in- 
teractions within the damage zone and proposes a link 
between the roughness exponent and the critical expo- 
nent v for the correlation length of the damage clusters. 
More generally, front propagation in random media has 
become a challenging problem related to the dynamics of 
interfaces in many different physic al sy stems theoretically 
connected, such as crack fronts [ill , m agnetic domain 
walls 01 > or wetting contact lines 0, Q3, where 
elasticity and disorder compete to shape the interface. 

In order to shed some light on the interactions between 
the crack front and material heterogeneities, a simplifi- 
cation to a two dimensional configuration -an interfacial 
crack- has been proposed both experimentally [Til Il7j 
and theoretically jlH l2Sl j . The interfacial configuration 
provides a higher resolution since all locations of the 
crack front belong to the same plane. Moreover, using a 
transparent material and a high resolution fast camera, 
the detailed complex crack dynamics can be captured, 
following the crack front with a high precision both in 



time and space |l9j . So far experiments have been fo- 
cused on the fracture front line morphology leadin g to 
the estimated roughness exponent £ = 0.55 ± 0.03 [la ], 
followed up by a longer study showing ( = 0.63 ± 0.03 
|l7j . First attempts have been recently performed to an- 
alyze the interfacial crack front dynamics [T^ . I20I ] . These 
studies have shown that the fracture front propagation is 
intermittent and can be described in terms of a Family- 
Vicsek scaling j2]| with a roughness C = 0.6 and a dy- 
namic exponent k = 1.2 ± 0.2. 

In this Letter, we study a system first studied exper- 
imentally by Schmittbuhl and Mal0y 0, ^| . Whereas 
previous studies focused on the morphology of the in- 
terfacial crack 0, we focus on the local crack dynam- 
ics, and on the distribution in both time and space of 
the waiting time during pinning events. To address this 
problem, we introduce a new analysis procedure in or- 
der to study the local waiting time fluctations. The im- 
proved experimental techniques and resolution allow to 
show that the dynamics of the fracture front is driven by 
local irregular avalanches with very large size and veloc- 
ity fluctuations, and anisotropic shapes whose scaling is 
directly linked to the self-affine scaling of the crack front 
itself. This new set of experiments also confirms earlier 
results on such systems 0, ^| . 

We describe here experiments where two Plexiglas 
plates are annealed together to create a single block with 
a weak interface ^(|. The plates are of dimensions: 
32cm x 14cm x 1cm and 34cm x 12cm x 0.4cm, and 
annealed together at 205°C under several bars of nor- 
mal pressure. Before annealing, both plates are sand- 
blasted on one side with 50/zm steel particles or 100/im 
glass beads. Sand-blasting introduces a random topog- 
raphy which induces local toughness fluctuations during 
the annealing procedure. We have measured the profile 
of a sand-blasted Plexiglas surface, using a white light in- 
terferometry technique (performed at SINTEF Oslo lab- 
oratory) and found that the local irregularities have a 
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FIG. 1: a) Typical example of a picture recorded by the high 
speed camera (Photron Ultima) during an experiment with an 
average crack front speed (v) = 28.1/wn.s -1 , and a pixel size 
a — 3.5/xm. The solid line represents the interface separating 
the uncracked (in black) and cracked parts extracted after im- 
age analysis, b) Gray scale map of the waiting time matrix 
deduced from 10000 front positions recorded at a rate of 1000 
fps. The darker parts, the longer waiting times, c) Spatial 
distribution of clusters (in white) corresponding to velocities 
ten times larger than the average crack front speed. 



characteristic size about 15/xm |22|. While the upper 
Plexiglas plate is clamped to a stiff aluminium frame, 
a press applies a normal displacement to the lower one 
(lcm thick) at a constant low speed which results in a sta- 
ble crack propagation in mode I [lr| . The fracture front 
is observed with a high resolution fast camera mounted 
on a microscope. Two different cameras have been used, 
a Kodak Motion Korder Analyzer CCD camera which 
records up to 500 frames per second (fps) with a 512x240 
pixel resolution, and lately a much more powerful one, 
a Photron Ultima CMOS camera. Using this camera at 
a spatial resolution of 1024 x 512 pixels, and an acquisi- 
tion rate of 1000 fps we can follow the stable crack front 
during more than 12s (recording up to 12288 images). 
Different experiments have been performed varying the 
acquisition hardware, the microscope magnification cor- 
responding to a pixel size between 1.7 to 10/im, and the 
average front line speed ranging from 0.35 to AOfim/s. It 
is important to note that in all cases, the pixel size is 
smaller than the size of the local irregularities of about 
20fim due to the sand-blasting process. 

In order to analyze the local waiting time fluctuations 
and the burst dynamics, we propose the following proce- 
dure: the fracture front lines extracted from image analy- 
sis of the digital pictures (see Fig.QJi) are added to obtain 
a waiting time matrix W(x, y). This matrix has the di- 
mension of the original image and an initial value equal 
to zero. We add the value 1 to the matrix element w 
corresponding to each pixel of the detected front line po- 
sition (x,y). This procedure is performed for all frames 



FIG. 2: The velocity distribution P(v/{v)) as a function of 
the scaled velocity v/{v) for different experimental conditions 
(various average crack front speeds (v) and pixel sizes a). A 
fit (solid line) to all data for v > (v) has a slope —2.55. Inset 
shows the corresponding waiting time distribution P(w/{w)) 
as function of the scaled waiting time w/{w). The solid line 
represents a fit to all the data for w < (w) with a slope 0.55. 



of a given experiment in order to obtain the final waiting 
time matrix W(x, y). A gray scale map of this matrix is 
shown in Fig. ^d. The spatially random toughness along 
the weak interface generates a rough crack line in pin- 
ning the crack front (Fig. and triggers a rich burst 
activity on a wide range of length scales. The numerous 
and various regions of gray levels suggest this intermit- 
tent dynamics (Fig-QJ)). It is important to mention that 
the image recording is so fast that there are basically no 
holes in the waiting time matrix W(x, y), i.e. no regions 
of zero values (apart from below the first front, above 
the last one, and a few artifacts due to impurities in the 
sample). Then, we can deduce from W(x, y), a matrix 
V(x, y) of the local normal speed of the interface at the 
time when the front went through a particular position, 
by computing the inverse value of the corresponding ma- 
trix element w of W(x, y) multiplied by the ratio of the 
pixel size a and the typical time between two images St. 
Therefore, we can associate to each pixel corresponding 
to the crack line in each image, a local front velocity 
v = — jz . Finally, we can obtain the probability distribu- 
tion functions of the local waiting time w and the local 
front velocity v, by estimating the occurrence number of 
each measured waiting time or velocity on all the pixels 
in all the fracture front line images. The velocity dis- 
tribution P(v/(v)) is shown in Fig. in a log- log scale. 
A data collapse is obtained for all different experimen- 
tal conditions by scaling the local velocity v with the 
average crack front speed (v) which varies from one ex- 
periment to another. A clear power law behavior of the 
velocity distribution P(v/(v)) ~ (w/(v)) _r ' is observed 
for velocities larger than (v) with a crossover to a slowly 
increasing function for velocities smaller than (v) . A lin- 
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ear fit to the experimental data for v/(v) > 1 gives a slope 
—j] = —2.55 ± 0.15. The inserted figure shows a double 
logarithmic plot of the corresponding waiting time dis- 
tribution P(w/{w)), where w is the waiting time, and 
(w) the average waiting time for each experiment. A lin- 
ear fit to the experimental data for w < (w) has a slope 
0.55 ± 0.15 consistent with the exponent r\ — 2 deduced 
from the velocity distribution. The power law distribu- 
tion of the local velocities confirms once again previous 
observations, revealing a non trivial underlying dynamics 
as observed on a fast video recording. It is important to 
note that even though the first moment of the velocity 
fluctuations (v) exists, the second and higher moments 
are ill defined and dominated by the largest velocity fluc- 
tuations. In an earlier work, the velocity distribution was 
investigated with a different method, based on the dis- 
tance between subsequent fronts at a given time interval 
|l9l |. However, such method proved out to produce re- 
sults depending on the time interval chosen. Indeed, a 
short time between the front only gave contribution from 
the high velocity part of the distribution while a long 
time between the fronts gave a peak around the average 
velocity only. By using the concept of waiting time intro- 
duced in the present Letter, we are able to measure both 
high and low velocities. In the present case there exists 
a well defined length scale a at which the velocity can 
be measured. Using different magnification of the micro- 
scope, we have checked the robustness of our procedure 
and shown the reproducibility of our results for different 
pixel sizes. 

In order to analyze the local burst activity, we consider 
a thrcsholdcd matrix generated from the velocity matrix 
V(x, y), by setting the matrix elements v equal to one 
for v > C ■ (v) and zero elsewhere, where C is a constant 
of the order of a few unities. Fig. QJ; shows the spatial 
distribution of clusters of different sizes obtained from a 
thresholded matrix with a threshold level C — 10. The 
white clusters correspond to velocities 10 times larger 
than the average crack front speed (v) . Then, we can ex- 
tract from this thresholded velocity matrix the size dis- 
tribution of the high velocity bursts. The clusters con- 
nected to the first and last front, and thus belonging to 
the upper and lower white parts are excluded from the 
analysis. 

On Fig. |21 the cluster size distribution P(S/(S)) is 
shown for different experiments at a given threshold value 
C = 3. First, we show a data collapse for all the dif- 
ferent experiments performed by rescaling the clusters 
size S with the average burst size (S). Moreover, we 
clearly observe that the burst size distribution P(S/(S)) 
follows a power law with an exponent 7 = 1.7 proving 
that the burst dynamics occurs on all length scales. We 
have checked that this critical behavior, and in particular 
the exponent 7 = 1.7 ± 0.1, is really robust: normalizing 
by the average burst size (S), we can rescale all the dif- 
ferent distributions corresponding to diverse experimen- 
tal conditions and a wide range of threshold level values 
2 < C < 20 (sec inset of Fig.©. 
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FIG. 3: Burst size S distribution P(S/{S)), normalized by 
the average burst size (S), for different experimental condi- 
tions (the various symbols correspond to those on Fia. ufy . The 
bursts detected for each experiment correspond to clusters of 
velocities 3 times larger than the average crack front speed. A 
fit on all the data (dashed line) gives a slope equal to 1.71. 
Inset: Normalized bursts size distribution P(S/(S)) averaged 
over all the different experimental conditions, for a wide range 
of different threshold levels C.A fit to all the data, cutting the 
largest clusters at which a cut-off appears due to the lack of 
statistics (solid line), gives a slope equal to 1.67. 



We expect a connection between the spatial scaling of 
the bursts, and the self-affine scaling of the front line 
itself on large scales. To investigate this, we have for 
each cluster S chosen the smallest rectangular bounding 
box enclosing it. The size of the bounding box gives the 
length scale L y of the clusters along the growth direction 
and the length scale L x of the clusters along the average 
front line orientation. Fig. 0] shows the dependence of 
the average size (L y ) on the length scale (L x ) in a dou- 
ble logarithmic plot, for a wide range of threshold values 
2 < C < 20, averaged over all the different experimen- 
tal conditions (different pixels sizes and average crack 
front speeds). We clearly see that the avalanche clusters 
become anisotropic above a characteristic length scale 
Ld ~ 15/im. This typical size corresponds to the cor- 
relation length for the disorder introduced by the sand- 
blasting technique j^. Below Ld the local toughness is 
marked by the same individual asperity and as a result 
the thresholded velocity bursts appear isotropic. A fit to 
the data points for L x > 15/j.m gives a slope H = 0.66 
consistent with previous independent estimates of the 
roughness exponent ( = 0.63 ± 0.03 for the fracture front 
line. This result shows that the system exhibits sclf-affinc 
scaling with the same roughness exponent £ for the local 
burst as the fracture front line and brings a new confirma- 
tion of the roughness exponent found in our experiment, 
which is different and higher than most present theoret- 
ical or numerical predictions [Tsl . I231 0> HH • 

As mentioned in introduction, the scaling behavior of 
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FIG. 4: Average length scale (L y ) as function of the length 
scale (L x ), for different threshold levels C, averaged over all 
the different experiments performed. The solid line is a fit to 
the data points for L x > 15/itm and has a slope H = 0.66, 
consistent with the roughness exponent of the fracture front 
line. The dotted line represents the curve (L y ) = (L a ) and 
serves as a guide for the eye. 

elastic interfaces in random media is involved in vari- 
ous physical systems. Actually, the discrepancy between 
theoretical and measured roughness exponent has also 
been reported recently, for contact lines of helium-4 p| 
and water meniscus |l4j ITEj propagating on rough sub- 
strates. It confirmed that the Joanny-De Gennes model 
|26[ usually proposed to describe the contact lines dy- 
namics, which leads to the same kind of e quat ion of mo- 
tion for crack fronts [53 i is n °t sufficient |l5l |. Interest- 
ingly the roughness exponent found respectively equal to 



0.56 ± 0.03 and 0.52 ± 0.04 are close to the fracture front 
line roughness. Besides, for the helium-4 meniscus, power 
law avalanche size distributions have been measured with 
exponents 7 from 0.99 to 1.3, depending on the contact 
angle, which is different from the exponent 7 = 1.7 found 
in our experiments. However, avalanches are defined by 
Prevost et al ^| based on subtraction between fronts, 
which is different from the present technique, and this 7 
exponent could be sensitive to such definition. 

Recent simulations based on a quasi-static model and 
interpreted as a stress weighted percolation problem [2^ 
give for the first time consistent results with the exper- 
imental roughness and dynamic exponent [l9l |20| . The 
insensitivity of the velocity distribution in our experi- 
ments on the average velocity of the front gives support 
to the quasi-static assumptions used in these simulations. 
However, the simulated process zone was not observed in 
our experiments above the micrometer scale. It should 
be mentionncd that d yna mical effects have also been in- 
troduced in a model [25j with a full elastodynamic de- 
scription where elastic waves may trigger instabilities and 
modify the roughness of the crack front, leading to the 
value C = 0.5 j2|. 

No theory or simulations so far have investigated the 
local velocity or the burst fluctuations. It will be of great 
interest to perform these analysis on the numerical mod- 
els for a direct comparison with our experimental work. 
This will hopefully clarify the importance of dynamical 
effects in modeling the fracture front propagation. 
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